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: ABSTRACT 

(N ■ 

Rotation and activity are key parameters in stellar evolution and can be used 
O ■ to probe basic stellar physics. Here we present a study of rotation (measured 

as projected rotational velocity v sini) and chromospheric activity (measured as 
Ha equivalent width) based on an extensive set of high-resolution optical spectra 
j^j ■ obtained with the MIKE instrument on the 6.5 m Magellan Clay telescope. Our 

targets are 74 F-M dwarfs in the young stellar associations rj Chamaeleontis, TW 
Hydrae, (3 Pictoris, and Tucana-Horologium, spanning ages from 6 to 30 Myr. 
While the Ha equivalent widths for most F and G stars are consistent with pure 
photospheric absorption, most K and M stars show measurable chromospheric 
emission. By comparing Ha equivalent width in our sample to results in the 
literature, we see a clear evolutionary sequence: Chromospheric activity declines 
steadily from the T Tauri phase to the main sequence. Using activity 
indicator, we find a plausible age range for the Tuc-Hor association of 10-40 Myr. 
Between 5 and 30 Myr, we do not see evidence for rotational braking in the total 
sample, thus angular momentum is conserved, in contrast to younger stars. This 
difference indicates a change in the rotational regulation at ~5-10 Myr, possibly 
because disk braking cannot operate longer than typical disk lifetimes, allowing 
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the objects to spin up. On timescales of ~ 100 Myr there is some evidence 
for weak rotational braking, possibly due to stellar winds. The rotation-activity 
relation is flat in our sample; in contrast to main-sequence stars, there is no linear 
correlation for slow rotators. We argue that this is because young stars generate 
their magnetic fields in a fundamentally different way from main-sequence stars, 
and not just the result of a saturated solar-type dynamo. By comparing our 
rotational velocities with published rotation periods for a subset of stars, we 
determine ages of 13^g Myr and 9l® Myr for the rj Cha and TWA associations, 
respectively, consistent with prev ious estimates . Thu s we conclude that stellar 



radii from evolutionary models by lBaraffe et al.l (119981 ) are in agreement with the 
observed radii within ±15%. 

Subject headings: stars: pre-main-sequence, chromospheres, evolution, magnetic 
fields, rotation 



Introduction 



About a decade ago, only a few nearby low-mass stars with ages between 5 and 30 Myr 
were known. Therefore, our knowledge about stellar evolution in this pre-main-sequence (and 
post-T Tauri) phase depended largely on interpolating between the well-studied clusters at 
ages < 5 Myr, e.g. the Orion nebula cluster (ONC), IC348, and NGC2264, and benchmark 
zero age main sequence (ZAMS) clusters with ages between 30 and 100 Myr, e.g. the Pleiades, 
IC2391, and a Per. In the last ten years, however, several nearby and young stellar associa- 
tions have been discovered, which provide us with the target samples for in-depth studies of 
the stellar properties in this critical age range. Most objects in these associations are spread 
over large areas of the sky, and have been identifi ed primarily based on satellite all-sky sur- 



vey data from ROSAT, IRAS, or Hipparcos fe.g. iKastner et al 



1997 



Mamaiek et al. 



Barrado y Navascues et al.lll999l ; IZuckerman fc Webbll2000l ; iTorres et al 
200lh. 



2000 



1999 



Zuckerman et al 



This paper concentrates on the r\ Chamaeleontis cluster (hereafter r\ Cha, age 6 Myr), 
the TW Hydra (TWA, 8 Myr), (3 Pictoris Mo ving Group (BPMG 12 M yr), and Tucana- 
Horologium (TH, ~ 30 Myr) associations (see IZuckerman fc Song! 120041 . for a review). As 
these associations are all located within 100 pc of the sun, and for all of them, a significant 
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population of low-mass stars has been identified within the past few years, they are ideal 
for probing stellar evolution from the T Tauri phase to the main sequence. Two stellar 
properties that are believed to undergo significant changes in this age range are magnetic 
activity and rotation, which are the main focus of this paper. 

Magnetic activity is a collective term for all phenomena caused by the operation of the 
stellar dynamo, including photospheric spots, chromospheric flares and plages, and coronal 
emission. It is observed mainly in the X-rays in the case of coronal activity and in op- 
tical emission lines for chromospheric processes. According to our current understanding, 
main-sequence, sun-like stars generate large-scale magnetic fields by means of an aQ-type 
dyn amo, which operates in t he transition layer between convective zone and radiative core 



(see iSchrijver Zwaanll2000l . and references therein). This type of dynamo is driven by the 
rotation of the star, and hence a strong rotati on-activity correla tion is observed for solar-type 
stars with ages >100 Myr. initially found bvlSkumanichl (119721) and subsequently con firmed 
by numerous studies (e.g. Patten fc Simon 19961 ; Terndrup et al. 2000 ; Barnes 2001 ). The 
correlation, observed both with coronal and chromospheric activity indicators, is usually 
seen as a linear increase of activity indicators with increasing rotation rate, followed by a 
saturation of the activity at high rotation rates. Young T Tauri stars, however, are fully con- 
vective, and thus cannot harbor a solar-type dynamo. As they evolve to the main-sequence, 
a change in the magnetic field generation and, as a consequence, in the magnetic activity 
properties and their connection to stellar rotation can be expected. 

The rotation of stars themselves changes critically in the pre-main-sequence phase. As 
the stars contract hydrostatically to the ZAMS, their rotation rates increase as a consequence 
of angular momentum conservation. In the first few Myrs of their evolution, the stars are 
believed to lose significant a ngular momentum due to magnetic interaction with their cir- 



cumstellar environments (see lHerbst et al.l 120071 . for a review), either by magnetic coupling 



betwe e n star and disk - a pro cess often referred to as 'disk locking', see for example iKoenigl 
(119911) ; Edwards et al.l (119931 ) - or by an accretion powered stellar wind (IMatt fc Pudritz 
20051 ). Once the disk is gone and the accretion has stopped, angular momentum removal is 
mainly controlled by stellar winds generated by magnetic activity. On the main-sequence, 
rotational braking by winds leads to a decline of both rotation and potentially (via a rota- 
tionally driven dynamo, see above) magnetic activity (jSkumanichl Il972l ; iBarnesI 120011 ). In 
summary, observational studies of magnetic activity, rotation, and their interrelationship can 
probe fundamental changes in stellar physics occurring between ~ 1 and 100 Myr. 

In this paper we present a study of rotation and activity in young stars with ages 
between 6 and 30 Myr, based on an extensive set of high-resolution spectra. By comparing 
the results from the associations listed above with literature data for younger and older 
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clusters, we investigate the pre-main-sequence evolution of these properties. For the first 
time, we analyze a comprehensive set of Ha line measurements together with spectroscopic 
rotational velocities v sin i for a large sample of objects in our fou r targe t regions. In some 
sense, this is complementary to the work by Ide la Reza fc Pinzonl (120041 ). who published a 
rotation/activity study for TWA, BPMG, and TH based on vsmi and X-ray data. 

In a previous paper, w e investigated disk accretio n for these associations based on the 
same multi-epoch dataset (jJayawardhana et al.ll2006l ). Our cleaned target sample for the 
current paper is introduced in £J2J In §31 we present our observations, the data reduction 
and spectral analysis, as well as a detailed assessment of the measurement uncertainties. 
Subsequently, we investigate Ha emission and rotational velocities (§5]) as a function of 
age and mass. In §5.31 we focus on the rotation/activity connection. By comparing vsmi 
with previously measured rotation periods, we obtain constraints on stellar radii for a subset 
of the targets (§SJ). We summarize our results in 



2. Target selection and properties 



This study is based on multi-epoch spectra of about 100 likely members of the four 
as sociations listed in jfi] From the total sample, we excluded objects identified as accretors 



m 



Jayawardhana et al. 



(2006), because accretion aff ects emission lines and thu s makes it 



difficult to assess chromospheric activity reliably. In iJayawardhana et al.l (120061 ) , the main 
criterion to distinguish between accretors and non-accretors is the linewidth of the Ha emis- 
sion feature. As thresholds we adopted 10 A for the equivalent width (EW) and 200 km s -1 
for the 10% width. We used the shape of Ha and additional emission lines like He 6678 A as 
complementary probes for accretion. Four objects (two in r\ Cha, two in TWA) are classified 
as accretors because their linewidths exceed both thresholds and they show He 6678 A in 
emission. One additional object, rj Cha 11, has very broad Ha as well. Although the EW is 
below 10 A, we classified it as an accretor based on a strong red-shifted absorption feature 
in Ha, a clear indication of infalling material. In total, five accreting objects were excluded 
from the rotation/activity analysis in this paper. 

Similarly, we excluded known unresolved binaries, since binarity might i ntroduce addi- 
tional line broadening. Many known binaries in these associations are listed by lZuckerman fc Song 
( 120041 ). a few more have been ide ntified in separate papers, e.g. the close AO-resolved binary 
TWA5A farandeker et al.ll2006f ). Additionally, the A-stars pPic and TWA11A as well as 
the subgiant HD 1555A were removed from the object list. The final cleaned target list 
comprises 74 stars. 
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The spectral types fo r our s ources have been collected from IZuckerman fc Song (120041 ) 



and Ide la Reza fc Pinzonl (120041 ) . In a few cases where spectral types are missing in the 
literature, we obtained an approximate spectral type by simply comparing the appearance 
of the spectrum with other stars from our sample. For plotting purposes, we convert spectral 
types to a linear numerical scale, where one unit corresponds to one subclass, and zero is 
assigned to O0 stars. Spectral type F0 then corresponds to 30 while M0 corresponds to 60. 

Our targets span the spectral type range from F5 to M3. The distribution of spectral 
types favours late-type objects, as the sample includes 36 M- and 19 K-stars. Spectral types 
are not equally distributed in all four associations: With one exception, all objects in rj Cha 
and TWA are K/M-stars. In the two older groups BPMG and TH, however, the distribution 
is balanced with more F/G-type objects. Effective temperatures have been derived by fitting 
models to the observed spectra, and the detailed results from this analysis will be published 
in a forthcoming paper (Mentuch et al., in prep.). In this paper, we make use of these T e g 
to determine Lna/^boi in §5.31 



The magnetic field generation and, as a consequence, magnetic activity and possibly 
rotation depend critically on the stellar internal structure. For pre-main-sequence stars, 
it is particularly important to clarify if they are still fully convective or if they already 
have developed to a solar-like structure with a radiative core and a convective envelope. As 
outlined in Jfl a radiative core is mandatory for the operation of a transition zone, solar-type 
dynamo, while in fully convective objects alternative dynamo mechanisms are at work. The 
age at which a radiative core appears is mass-dependent: While solar-mass stars are only 
fully convective for 1-3 Myr, stars with 0.5 M Q need 10-20 Myr to develop a radiative core. 
Objects with masses below 0.35 M® are fully convective throughout their evolution. 



We used the theoretical evolutionary tracks by lChabrier fc Baraffd (119971 ) and lD'Antona fc Mazzitelli 



(11994 updated version 1998) to assess the internal structure of our targets. In the three 
youngest associations, rj Cha, TWA, and (3 Pic, all objects with masses > 0.6 M are not 
anymore fully convective, according to the models, where the depth of the convection zone 
increases towards higher masses. For TH, the mass limit between radiative and fully convec- 
tive is roughly at 0.3-0.4 M , i.e. comparable to the main sequence. Converted to spectral 
types, this implies that all objects in our sample with spectral types earlier than M0-M2 
already have a radiative core. In TH, this applies to all objects earlier than M3. All these 
limits are subject to uncertainties due to model inconsistencies, age spread, and uncertainties 
in the mass/temperature/spectral type conversion. 

This essentially means that almost all our targets in i] Cha and TWA are still fully 
convective or just at the limit to develop a (small) radiative core, while in (3 Pic, six out of 
16 objects are clearly already in their radiative evolution. In TH, with its large fraction of 
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K/G/F stars, almost all objects have a substantial radiative core according to the models 
(i.e. basically a solar-like internal structure); only the three M3 stars may still be fully 
convective. We will use these considerations as a guideline when discussing evolution in 
magnetic activity and rotation. 



3. Observations, spectral analysis, uncertainties 

3.1. Observation and Data Reduction 

The observations were taken on 12 nights distributed over four observing runs between 
December 2004 and July 2005, using the echelle spectrograph MIKE at the Magellan Clay 
6.5 m telescope on Las Campanas, Chile. In total, the dataset comprises of ~ 650 high- 
resolution spectra with each star from our sample being observed on average over 6 times. 
The availability of multi-epoch data for each star allows for a reliable assessment of variability 
in the emission lines and provides us with high red undancy in the v sin i measu rements. For 



more details about the observing runs, see the log in lJayawardhana et al.l (120061 ). In addition 



to our target stars, we observed a sample of ~ 2 slowly rotat i ng sta ndard stars covering 



the spectral range from F8 to M2 from the list of iNidever et al.l (120021 ) . 



MIKE is a double echelle slit spectrograph, consisting of a blue and a red arm. In this 
study, we only make use of the red part, with spectral coverage from 4 900 to 9 300 A. With 
a / . / 35 slit width and no binning, our spectra have a resolution of R ~60 000. Depending 
on the brightness of the target, the integration time was up to 30min. To accommodate for 
the slanted spectra and the wavelength dependence of this tilt, we developed a customized 
software package in ESO-MIDAS for data reduction. The details of this reduction procedure 
will be discussed in a forthcoming paper (Brandeker et al., in prep.). 



3.2. Spectral Analysis 

As an indicator of chromospheric magnetic activity, we measured the EW of the Ha 
feature at 6 562.8 A in the multi-epoch spectra in our sample. Depending on the spectral 
type and the level of chromospheric activity, Ha is seen either as emission or absorption in 
our spectra. The EWs are obtained by integrating the Ha line after continuum subtraction; 
absorption is defined as positive EW, emission as negative. The continuum is approximated 
by a linear fit to data points in small regions to the immediate right and left of the feature. 
These continuum defining regions are between 6 A and 10 A in width, which corresponds to 
200 to 450 data points in the spectra. 



- 7- 



We determined the projected rotational velocity vsini by creating a template from 
spinning up a slowly rotating standard star of similar spectral type to the target star. The 
differences in spectral types between target and standard were usually less than two sub- 
classes. The templates are created in lkms -1 increments; we adopted as vsini the value 
which produces the best x 2 fit between target and template. This procedure was carried out 
in three spectral regions for every available spectrum of a given star, with the average value 
taken to be the final vsini of the target star. The three spectral regions used for the fit 
have a width of 30 A and are centered at 5310, 6250, and 6815 A. They are chosen to avoid 
emission features and transition regions between different orders of the spectra. Since our 
earliest spectral type standard is an F7 star, we are unable to obtain v sin i measurements for 
stars that have spectral types significantly earlier than this. The availability of multi-epoch 
spectra allows us several independent measurements per star for both Ha EW and vsini. 
The standard deviation from the average provides an estimate of the consistency of our re- 
sults and the variability in Ha. All average Ha EW and vsini measurements are listed in 
Tables [M 

For a subset of our objec ts, rotational velocities hav e already been published in the 



literature, compiled recently by lde la Reza fc Pinzonl (120041 ) . We use these values as a control 



sample in the error analysis (see §3.3p . 



3.3. Error Analysis 

Strong emphasis was placed on the reliability of the rotational velocities. Particularly, 
we are interested to find the lower limit in v sin i for which our results are trustworthy. 
We analyze the reliability and the accuracy of our measurements in a variety of ways. As 
explained in §3.21 we measure vsini in three different wavelength bins in each available 
spectrum of our multi-epoch dataset, giving us at least 3 and on average 15 independent 
measurements per target. The value adopted for the star is the average over these individual 
measurements. From the individual measurements for each star we also obtain the standard 
deviation (of the average EW), a measure for the spread in our datapoints for a given star, 
and thus a probe for the internal consistency of the method. Fig. [1] shows the standard 
deviation a versus the average value of vsini for all objects. As can be seen in the figure, 
a is less than 2.5kms _1 for the majority of our objects, confirming the reliability of our 
method. Only 4 out of 74 stars, or 5%, have standard deviations exceeding 5kms _1 . Thus, 
a conservative estimate of the internal accuracy in our vsini values is 5kms _1 . We adopt 
this as a preliminary lower vsini limit. The standard deviation a is used in the remainder 
in the paper as the uncertainty in individual v sin i measurements. 



- 8- 



In the next step, we compare the measured vsini with published values, if available. As 
shown in Fig. [2j the majority of our values are within 8kms _1 of those obtained from the 
literature. The two stars for which the difference between the measured and literature v sin i 
measurements exceed this amount are HIP 2729 and PZ Tel, which are both fast rotators 
with vsini exceeding 60kms _1 . Although the absolute difference between literature and our 
v sin % for both of these stars is large, the relative uncertainty is comparable to the slower 
rotating objects in the sample. Thus, the comparison with the literature confirms the validity 
of our measurements. Since this comparison encompasses both the errors in ou r own v alues 



and the errors in the literature values (which are at least 3kms _1 , see iReidl (120031 )). we 
find that the absolute average uncertainty in our vsini measurements is comparable to the 
internal accuracy, which has been found to be 5kms _1 (see above). Thus, we can safely 
trust all values down to at least 5kms _1 . 

Finally, we explore the reliability of the method by using a standard star (called A in 
the following) to measure the vsini of another standard (called B) of similar spectral type 
that has been spun-up to a determined value. That way, we produce a template with known 
rotational velocity, for which we can test the accuracy. This assumes that our standards are 
slow rotators with rotational velocities well below our detection limit. Standard B is spun-up 
in steps of lkms" 1 between and lOkms -1 . We repeat the measurements reversing the 
roles of the standards, i.e. switching A and B. These sets of measurements are done for four 
pairs, representing cases in which we have more than one standard star per spectral type. 
This produces eight different pairings of standards with the same spectral type. In six of 
these pairings, the measured v sin i approaches zero when the imposed v sin i approaches 
zero, demonstrating that in these cases our method produces highly accurate measurements 
of v sin In two pairings, however, the measured vsini levels off at 5kms _1 for imposed 
vsini < 3kms _1 . The best explanation for this result is that in these cases the standard 
A has a vsini that is 5kms _1 greater than that of standard B, and thus we are never 
able to measure values lower than 5kms _1 . This is evidence for non-negligible rotational 
velocities in some of our standards, which is probably the dominant source of uncertainty in 
our measurements. The test shows, however, that these uncertainties are unlikely to exceed 
5kms _1 , thus confirming the lower reliability limit for our method. In the following, we 
treat all vsini values lower than 5kms _1 as upper limits. We note that this affects only five 
objects. 

We estimate the uncertainty in the measurements of Ha equivalent width to be on 
average 0.2 A. This estimate was arrived at by determining the average value of equivalent 
width obtained from regions of the continuum close to Ha that do not contain visible emission 
or absorption features. As can be seen in Tables [THU the standard deviation a in the average 
Ha EW, calculated over the multi-epoch data, is in many cases clearly higher than the 
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measurement error, indicating significant variability (see §4.2j) . 



4. Magnetic activity 

Magnetic activity and its interplay with rotation is a complex problem, as it depends 
on stellar age, mass, interior structure, and possibly interactions with disks in the early 
evolution. To disentangle the different processes, we will start in this section by analyzing 
Ha EW as a function of age and mass, We note that many of our targets have been identified 
in X-ray surveys, thus the selection might be biased towards more active objects. Therefore, 
in the discussion of activity, we prefer to use criteria based on the upper limit of activity in 
our sample rather than the lower limit, since the latter one might be biased. 



4.1. Evolution of chromospheric activity 

The Ha feature is used routinely as an indicator of chromospheric activity, originating 
from photoionization and collisions in the hot chromosphere. A main tool to investigate the 
chromospheric activity is to plot Ha EW vs. effective temperature, in our case represented 
by the spectral type. This plot is shown in Fig. [3] for our young targets. 

Two main features are obvious from this figure: a) There is apparently no difference in 
the distribution of datapoints for our four target regions, spanning an age range from 6 to 
30 Myr. This is further strengthened by the fact the the distribution of EWs are statistically 
indistinguishable in the four regions. Thus, the activity levels in the four associations are 
faily similar, as far as we can tell with our data (see below), although their ages are somewhat 
different, b) Ha EWs are a strong function of spectral type. While mid F-type stars exhibit 
Ha absorption of ~ 5 A, the feature switches to emission at K2- K4 spectral types. Around 



MO, corresponding to masses of 0.7-0.8 M (IBaraffe et al.lll998l ). there is a clear 'knee' in 



the distribution; at the same time, the spread in EW increases, and the stars reach emission 
levels between and —10 A at early/mid M spectral types. 

The strong change of Ha EW with spectral type does not only reflect a change in 
chromospheric activity, as the EW is additionally affected by the drop in the continuum 
level with stellar luminosity and the photospheric absorption in Ha, which is about zero for 
M dwarfs and increases towards earlier spectral types. The combined effects of photospheric 
continuum drop and Ha absorption are estimated from the Ha EW for our standard stars, 
whi ch we already used as rotational velocity templates and which are selected to be non-active 



; sec 



Nidever et al.l 120021 ) . A linear fit to their EW as a function of spectral type is shown in 
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Fig. [3] as a dotted line; the object-to-object scatter around this line is typically ±0.5 A. We 
note that this dashed line is consistent with the EWs of non-active stars in the Hyades a nd in 



the field for late K and early M spectral types (H erbst fc M iller 



with published EW for field F and G dwarfs (IPeat 



1964 



1989 



Reid et al 



Strauss k, Ducati 



i n line with theoretical predictions for Ha EW without chromosphere from 



1995|), and 



19811). It is also 



Cram fc Mullan 



( 119851 ) . Thus, this line is an estimate for the pure photospheric contribution to the Ha EW. 



As can be seen in the figure, the dotted line follows the lower envelope of the EW for 
our young target stars between F8 and M2. All EWs for spectral types earlier than G8 are 
in agreement with the dotted line and thus pure photospheric values; thus, these stars do not 
show measurable chromospheric activity. Starting at spectral type G8, the EWs measured 
for our targets show increasing excess with respect to th e photospheric values, i ndicating a 
contribution from magnetic activity. As pointed out by I Cram fc Mullan! (119791 ). the onset 
of chromospheric activity will first tend to deepen the absorption feature by as much as 
0.5-1 A in relatively cool and thin chromospheres, where line formation is dominated by 
photoionization and not by collisions. This effect is not seen in our data, no object shows 
significantly more absorption than the photospheric values. Thus, objects in transition 
between essentially non-active chromospheres to the collision-dominated regime are rare. 
The maximum level of magnetic activity and the fraction of active objects increases rapidly 
from early K to mid M spectral types. 

These results can be compared to studies of older and younger objects in a similar 
spectral range. We use three criteria: 

A) The spectral type (or color) at which Ha changes from absorption to emission. This 
value was introduced by lHawley et al.l (119991 ) as an indicator of stellar age, as it is steadily 
shifting to later spectral types as the objects get older. In our sample, the transition is at 
spectral types K2-K4, but it is only accurately defined for objects in TH. As summarized in 
Fig. 5 of lHawley et al.l (119991 ). the transition occurs at early M types in the Hyades (age 0.5- 
1 Gyr), at late K types in the Pleiades (age 125 Myr), and at mid K types in IC 2602/2391 
(age 30-40 Myr). For objects in the 1-5 Myr age ra nge (i.e. you nger than our sample), the 
transition occurs at spectral types earlier th an K0 (iDahml 120 05) - in fact, objects without 
Ha in emission are very rare at these ages (IPoncet et al.lll998l ). Thus, the stars in TH fit 
nicely in the evolutionary sequence defined in the literature, indicating a steady decline of 
activity in the pre-main-sequence evolution. Using this criterion as an age indicator, we find 
that ages in the TH assoc i ation are most likely b e tween 10 and 40 Myr, confirming previous 
estimates by lTorres et al.l (120001 ) and ISong et al.l (120041 ) . 



B) The fraction of K and M type objects with Ha emission. In our sample, practically 
all K/M objects are above the photospheric values, indicating activity, which is also the case 
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for (non-accreting) stars younger than 6 Myr. A close to 100% fraction of acti ye stars is also 



seen in pre- main-sequence clusters like IC 2602/2391 with ages of 30-40 Myr (jStauffer et al. 



1997bl ). In contras t, only a small f raction of active stars is seen in the Hy ades at an age of 
0.5-1 Gyr (< 30%. lReid et~allll995h and in the old field population (< 10%. lHerbst fc Miller 
19891 ). Criterion B thus confirms the drop in activity for objects older 30 Myr, but some of 
this effect might be due to the selection bias in our sample towards highly active objects (see 
above) . 

C) The maximum level of activity in M type objects. M stars are the most active objects 
in our sample with EW ranging from zero to 11-12 A. To avoid being biased too much by a 
few extremely active objects, we do not take into account the most active 10% of the objects 
and thus obtain ~ 9 A in our sample. We note that this value is mostly determined by 
objects in the three younger regi ons. For com parison, the same procedure gives ~ 9 A for 
objects younger than 6 Myr (e .g. lDahmll2005l ). ~ 6.5 A in the Pleiades, and ~ 6.0 A in the 
Hyades (ITerndrup et al.l 120001 ) . Again this criterion confirms a decline of activity, which is 
particularly significant between 30 and 100 Myr. 

From all three tests, a clear evolutionary sequence is apparent. Our targets are signifi- 
cantly more active than objects in the Pleiades, Hyades, and older samples. Thus, there is 
a clear decline of chromospheric activity from 30 Myr to 100 Myr and beyond. This decline 
leads to a reduction of the activity level in all spectral types, which includes the complete 
disappearance of measurable chromospheric emission at earlier spectral types. As the objects 
age from 30 Myr to 1 Gyr, the spectral range of non-active objects extends to later spectral 
types. Eventually, among evolved main-sequence stars only a fraction of M-type objects can 
maintain chromospheric activity. The best interpretation of this behavior is a mass depen- 
dence in the lifetime of active chromospheres: According to the numbers given above, the 
timescale on which the chromosphere provides significant Ha emission are < 10 Myr for F 
and G stars, 10-100 Myr for K stars, and > 500 Myr for M stars. A similar decline of activity 
can probably also be seen in the flare frequencies (see §4.2p 

Maybe the best approach to explain this result is based on the different internal structure 
of the stars in our sample: While the F and G stars already have substantial radiative cores, 
most M stars are still fully convective. It is conceivable that once a rotationally driven au 
dynamo is able to operate at the transition between convective and radiative zone, angular 
momentum loss through stellar winds quickly shuts down the magnetic activity. Thus, 
the connection between dynamo activity and rotation could potentially be the origin of the 
quick disappearance of the activity for early-type objects. An alternative explanation for the 
decline of chromospheric activity on the pre-main-sequence might be that the change in the 
interior structure alters the properties of the magnetic surface field. This in turn might lead 
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to a change in the dominating mechanism responsible for the heating of the chromosphere, 
which will then affect the Ha emission. 

As already mentioned above, we do not see significant differences between the four 
associations covered with our sample. However, the spectral type coverage in the four regions 
is different, and thus it is not possible to carry out a more rigorous comparison of the four 
regions using the three criteria defined above. In TH, for example, the number of M-type 
objects is very low, hampering a reliable assessment of the fraction of active objects and their 
maximum activity level. On the other hand, the younger groups 77 Cha, TWA, and BPMG 
lack objects to spectral types between G5 and K5, compromising an analysis for the onset 
of activity in those regions. Thus, from our data alone we cannot definitely rule out activity 
evolution between 6 and 30 Myr. When comparing our data with younger objects, however, 
we see evidence for activity evolution on this timescale, in the sense that the transition to 
emission occurs at somewhat earlier spectral types in T Tauri stars. Taken together, the 
analysis in this section indicates that chromospheric activity steadily declines as the stars 
evolve from the T Tauri phase to the main sequence. 



4.2. Variable chromospheric activity 

Since we have more than one epoch for most of our targets, we are able to probe 
variability in the Ha emission. Because both photospheric Ha absorption and bolometric 
luminosity are not expected to change significantly (i.e. more than a few percent) for these 
objects, variability in Ha EW basically traces changes in the level of chromospheric emission. 
For a few objects in the youngest regions, weak levels of episodic accretion cannot b e excluded 



and might contribute somewhat to the variability (see I Jayawardhana et al.ll2006l ). 



The primary estimate of variability is the standard deviations in our EW time series. In 
Fig. H] (left panel) we plot the absolute values of Ha EW a vs. spectral type. The dashed line 
marks the measurement uncertainty. As can be seen from this plot, many objects with late 
spectral types show significantly higher Ha variations than expected from the formal error, 
indicating variability in activity. Interestingly, the onset of measurable variability occurs at 
early K spectral types, where Ha changes from absorption to emission. This confirms that 
the variations can indeed be attributed to chromospheric activity - stars without measurable 
activity and thus only photospheric Ha do not show variability. The plot shows no significant 
difference between the four groups, indicating that the level of variability does not strongly 
change between 6 and 30 Myr. 



Ha emission originates from active regions in the chromosphere, which are typically not 
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uniformly distributed. Thus, one main cause of the Ha variations is rotational modulation. 
Additionally, the light curves can be affected by flare activity and overall changes in the 
activity level, e.g. due to an activity cycle. Our time sampling makes it difficult to distinguish 
between these three scenarios. In most cases, we have only one spectrum per night per target; 
the longest time baseline is eight months. Rotational changes occur on timescales of the 
rotation periods, which are typically a few days for our targets. These changes are periodic, 
but with our sparse sampling we are not able to recover the periods. General activity level 
changes are a long-term phenomenon, and thus might introduce a gradual trend in our time 
series. Isolated flare events would be detectable, but only if they are clearly stronger than 
all other sources of variability. If several flares are present in our time series, it would again 
be difficult to identify the source of variability. 

We checked all Ha EW time series for signs of isolated flare events. Since the typical 
flare length in the optical wavelength range is at most a few hours, a flare would appear as a 
single Ha measurement with significantly stronger emission than all other datapoints in this 
particular time series. As a clear flare event, we accept a positive 3a outlier in the time series. 
It turns out that none of our objects exhibits such an event, although about 10 of them show 
2a outliers (for example, the active stars TWA 10 and AUMic in BPMG). Thus, strong, 
isolated flares are rare in our sa mple. As already mentioned, flares last typically 1-2 h in the 



optical (IGuenther fc Ball! 119991 ) . Given our conservative detection limit, however, we would 
detect them only in the first 10-20 min, when their effect is most pronounced. In total, we 
have about 400 spectra, which thus cover about 100 h. Therefore, the flare frequency derived 
from our spectra is < 0.01 hT 1 . Assuming a flare duration of ~ 2h, this corresponds to a 
flare rate of < 2%. 

There are few rel i able s tatistical constraints on (average) flare frequencies in the opti- 



cal. IGuenther fc Ball! (119991 ) derive rates of chromospheric flares for non-accreting T Tauri 
stars (ages ~ 2 Myr) and ZAMS stars (ages ~ 50 Myr) based on multi-epoch multi-object 
spectroscopy, applying a similar criterion as we have used, but for H/3 instead of Ha. They 
find flare frequencies of 0.06 h~ x for T Tauri and 0.006 h~ x for ZAMS stars, concluding that 
the average flare frequency drops by a factor of ten as the stars evolve from 2 to 50 Myr. Our 
result of < 0.01 h~ x is clearly lower than the value derived for T Tauri stars, which might 
indicate that our targets are in an intermediate evolutionary stage between T Tauri phase 
and ZAMS in terms of their flare activity. 
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5. Stellar rotation 

Stellar rotation is known to be a function of mass, age, and magnetic activity. These 
dependences will be discussed separately in the following subsections, with the goal of dis- 
entangling the involved processes. 



5.1. Rotation vs. spectral type 



Rotation is known to change as a function of stellar mass, mainly because the efficiency 
of angular momentum removal depends on magnetic activity, which in turn depends, as 
discussed in §4.11 on stellar mass. In Fig. [5] we plot v sinz versus spectral type, which we use 
as an indicato r of stellar mass. E arly K spectral type roughly corresponds to 1M , early 
M to 0.5 M (IBaraffe et al.lll998l ). The majority of the objects have rotational velocities 
below 60kms _1 , the four exceptions, (5 ±3)% of our total sample, are PZ Tel and in BPMG 
and HIP 108422, HIP 2729, and CD-53544 in TH. Objects with v sin i > GOkms" 1 are called 
ultrafast rotators in the following. 

The overall appearance of this plot is comparable to v sin i distributions in young clus- 
ters. In the ONC, for example, typical values for v sin i for G-M spectral typ es are in the 
range of 12-30 km s -1 , while higher mass stars tend to rotate somewhat faster ( Rebull et al. 



20021). For F-M spectral types, velocities >60kms 1 are in general rare (5-10%. ISicilia-Aguilar et al 
20051 ). consistent with our dataset. The usini distributions in ZAMS clusters like the 



Pleiades ( iTerndrup et al.l 120001 ) or IC 2391/2602 shows the same phenomeno logical appear- 
ance. In these clusters, the number of ultrafast rotators might be somewhat higher (~ 15%), 
as expected as a consequence of pre-main-sequence contraction and thus rotational acceler- 
ation (see §5.21) . 

For the early-type stars in our sample both the upper and the lower limit of the v sin % 
distribution decline steadily with spectral type. Excluding the ultrafast rotators, the upper 
limit decreases from ~ 50kms _1 at F5 to ~ 15kms _1 at K5, while the lower limit drops 
from ~25kms~ 1 to the detec tion limit of 5 km s" 1 in the same spectr al range. A similar 
trend is seen in the Pleiades ( IQueloz et al.lll998t iTerndrup et al.ll2000l ). There are at least 
two possible explanations: 

a) The timescale on which the rotation of the stars is braked as a consequence of star- 
disk interaction (see §5.21) . depends on spectral type, in the sense that early type objects 
loose their disks faster than later types. Evidence for mass- dependent disk lifetimes has been 
found recently ([Carpenter, et. al.l 120061 : IScholz et al.l 120071 ). but further tests are needed to 
clarify the impact on rotational evolution. 
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b) The effect can also be understood as a consequence of a change in the stellar interior 
structure: As already discussed in £|2] and §4.11 all objects earlier than MO in our sample 
do have a radiative core and thus are able to operate a solar-type dynamo. For these 
objects, a deep convection zone enables efficient ang ular momentum removal due to stellar 
winds and/or disk- locking (jSchrijver fc Zwaanll2000l ). At any given age > 5 Myr, stars wit h 
spectral types K have deeper convection zones than F-G stars (ID'Antona fc MazziteUilll994l ). 
Thus, as we approach later spectral types and the convection zones in the stars become 
progressively deeper, the rotational braking becomes more effective, resulting in reduced 
rotational velocities, as seen in Fig. [5j 



5.2. Rotational evolution in the pre-main-sequence phase 

To examine the evolutionary effects more in detail, we plot v sin i vs. age in Fig. [6j In the 
upper panel, we show only the four associations. As can be seen in the plot, the upper limit in 
v sin i increases with age; the ultrafast rotators are only seen in older associations. This trend, 
however, might be a result of small number statistics. We compared the distributions of v sin i 
using a double-sided Kolmogoroff-Smirnoff test. Specifically, we tested the null hypothesis 
'the v sin i distribution in two associations is the same'. It was found that with two exceptions 
all possible combinations of r] Cha, TWA, BPMG, and TH give likelihoods for the validity 
of the null hypothesis larger than 25%. When comparing TWA with older associations there 
is some marginal evidence for statistical differences, with false alarm probabilities of 6.2% 
(BPMG) and 6.3% (TH). In general, however, th e four datasets are fairly similar. This is 



consistent with the results of IStassun et al.l (119991 ) , who find similar v sin % distributions for 



the ONC (lMyr) and the Pleiades (125 Myr). Thus, the overall distribution of rotational 
velocities does not appear to change significantly in the pre-main-sequence phase. Please 
note that this does not necessarily imply consistency with conservation of angular momentum 
throughout the pre-main-sequence phase, as the objects undergo a strong contraction (see 
below for a more detailed assessment). 

A large scatter of rotation rates is seen at all ages. While the projection factor sini, age 
spread, and v sinz uncertainties all contribute to the scatter, the major reason for the large 
spread of the distribution is probably the spread in the initial rotation periods. In cluster s 



with ages of 1-2 Myr, the periods range from fractions of a day to ~ 20 d (IHerbst et al.ll2007l ). 
corresponding to rotational velocities ranging from < 5kms _1 to > lOOkms -1 . Moreover, 
the v sin i distribution (as well as the distribution of log v sin i) is highly asymmetric, hamper- 
ing a rigorous statistical analysis. To mitigate this problem when investigating the rotational 
evolution we work in the following with typical (median) v sini values for a given age, rather 



-16- 



than with individual datapoints. Please note that by averaging over the rotational velocities 
in one particular group, we loose any information about the spectral type dependence of the 
rotation, which has been discussed in §5.11 

In Fig. [6j we overplot the median values for each association as large octagons. We will 
compare these median vsini with simple models for the rotational evolution. As starting 
value for the models, we used the typical vsini of 8-15kms _1 (average 11.5kms _1 ) at ~ 



5 Myr given by IRebull et al.1 (I2004T). To take i nto account the pre-main-sequence contraction, 
we use radii from IChabrier fc Baraffd (119971 ) for a stellar mass of 0.8 M Q , which is typical 
for our sample. 

In the upper panel of Fig. [61 we plot the expected rotational evolution for two extreme 
cases, constant angular momentum with a solid line (model A ) and constan t angu lar velocity 
with a dashed line (model B). In this approach, we follow iRebull et al.l (12004 ) who have 
done a similar comparison for stars with ages from 1-10 Myr. In model B, the period is 
constant, as expected in a scenario with ideal 'disk-locking', and thus vsini oc R. Model 
A, on the other hand, shows purely the spin-up due to contraction and thus vsini oc R -1 . 
While both models are in good agreement with observations until ages of ~ 10 Myr, only 
model A is clearly consistent with the median vsini at 30 Myr. Model B, however, gives 
too low values for ages >10Myr; it truncates the vsini distribution at the 20% quartile in 
BPMG and at the 10% quartile in TH. Thus, from 10 to 30 Myr the objects show rotation 
rates rather consistent with conservation of angular momentum than with constant rotation 
period. Thus, the dominating effect for the rotational evolution in this time window is spin- 
up due to the pre-main-sequence contraction. This result is robust against uncertainties in 
the stellar radii, because only the ratio of radii is used in the calculation. 

In strong contrast to our finding, for ages < 5 Myr the rotational ey o lution closely 
follows the track for constant angular velocity, as concluded by IRebull et al.l (120041 ) . There 
is growing evidence for a strong rotatio nal braking in the first few Myr, most likely produced 
by interaction with accretion disks (e.g. Herbst et al. 2002 ; Rebull et alJbood) and preventing 



the stars from spinning up by essentially lo cking the rotation period (e.g. IRebull et al.ll2002 
Tinker et al.l 120021 ; iHerbst fe Mundtl 120051 ) . Our results now demonstrate that while the 
period may be locked until ages of ~ 5 — 10 Myr, in the following ~ 20 Myr the stars spin 
up without clear evidence for rotational braking. Thus, rotational acceleration (measured 
in period) becomes significan t at ages of 5-10 M yr - which is consistent with the typical 
lifetime of circumstellar disks (IHaisch et al.ll200ll ). Specifically, it has been shown that many 
of the youngest stars in our sample (i n r\ Cha and TWA) are affected by inner disk clearing 
measured from mid- infrared excess (IHaisch et al.l 120051 ; iJayawardhana et al.lll999l ). while 
the oldest objects (in TH) do not show any evidence for disks at mid-infrared wavelengths 
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( IMamajek et al.ll2004j ). Thus, the change of the rotational regulation at 5-10 Myr coincides 
with the disappearance of the inner disks. It has to be emphasized, however, that all these 
considerations only apply to the typical evolution. For individual objects, the period-locking 
timescale can vary by a lot - possibly due to different disk lifetimes. 

To follow the evolution to the main sequence, we compared our dataset with the rota- 
tional velocity data in the Pleiades. In the lower panel of Fig. [6j we plot the median v sin % 
for F to M stars (large octago n) together wi t h the quartile values (horizontal bars). These 



numbers have been taken from lQueloz et al.l (119981 ) (their Fig. 6, averaged over all masses). 



In this plot we show for each model two evolutionary tracks, the first starts at 6 Myr and 
calculates forward in time (as in the upper panel), the second starts at 125 Myr and calcu- 
lates backwards. Solid lines show again model A, i.e. conservation of angular momentum 
without any rotational braking. The tracks from model A are barely consistent with the ob- 
servational data. When started at 5 Myr, the predicted median in the Pleiades is 21kms _1 
and thus too high; when started at 125 Myr, they give a median of 5kms _1 at 5 Myr, which 
is clearly too low. Thus, rotational braking is likely involved in the evolution to the ZAMS. 

On the main-sequence, rotation is mainly braked by angular momentum losses due 
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momentum losses according to the Skumanich law, again calculated in both directions. The 
tracks from model C, however, are clearly not in agreement with the observations. When 
calculated forward, the predicted median for the Pleiades is well below the detection limit. 
Conversely, for 5 Myr the model gives an unrealistically high median. Thus, Skumanich 
braking appears to be too strong. We can reproduce the wsini evolution either by using 
an exponent of —0.1 to —0.3 instead of —0.5 in the braking law, by using an exponential 
braking law with v sinz oc exp (— t) or by switching on the braking at about half way through 
the pre-main-sequence evolution. The latter scenario is not implausible, as most objects in 
the considered mass range develop a radiative core and thus the pre-requisite to operate a 
solar-type dynamo after about 30 Myr, see ^TJ 

Thus, our comparison with models gives the following results: a) On timescales of 
~ 100 Myr, weak rotational braking, possibly due to a Skumanich-type activity-rotation 
connection, is required to find a good match to the observations, b) From 5-30 Myr the 
rotational evolution is fully consistent with angular momentum conservation; effects of pos- 
sible rotational braking are too weak to affect the v sin % distribution significantly. Again it 
should be emphasised that these results do only apply to the total sample. In §5.11 we do 
find that rotational velocities depend on spectral type for objects earlier than M2. Thus, for 
objects with ages between 5 and 30 Myr, stellar mass is the major factor which determines 
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the rotation, rather than age. 

In the previous sections we have already made connections between rotation and activity, 
to explain the evolution and mass-dependence of Ha emission and rotational velocities. The 
obvious next step is to investigate directly possible correlations between rotation and activity, 
which is the focus of the next subsection. 



5.3. The rotation-activity connection 



In order to obtain a physically meaningful picture of a possible connection between ro- 
tation and activity, we derived relative Ha luminosities (i.e. Ln a /Lb \) from the measured 
Ha EW. We focused on the objects with clear chromospheric Ha emission, and therefore 
excluded stars with spectral type earlier than K2 (see §4.ip . In a first step, we corrected 
the EW for photospheric absorption, using the correlation between photospheric Ha absorp- 
tion and spectral types derived in §4.11 from non-active reference stars (see dotted line in 
Fig. [3]). Objects with corrected EW < 0.5 A and thus insignificant chromospheric emission 
were excluded. The continuum at the wavelength of Ha was estimated using the STAR- 
dustyl999 model sp ectra, which are ba sed on the NextGen models refreshed with new water 
and TiO opacities ( Allard et al.ll2000i ). We measured the continuum flux at 6562 A for ef- 
fective temperatures ranging from 3 000 to 5 000K and logg = 4.0 by approximating the 
spectrum around Ha with a linear fit. This value was divided by the bolometric luminosity 
for the respective effective temperature. As a result, we obtain scaling factors as a function 
of effective temperature to convert the Ha EW to L HQ /£boi- Pl ease n °te that this conversion 
depends neither on the radii of the objects nor on the distances, which are poorly constrained 
for many of our targets. The effective temperatures for our targets will be published in a 
forthcoming paper, see £J2J 

Fig. [7| shows the relative Ha luminosities as a function of v sin i. Please note that by 
excluding non-active (earlier type) objects, the clear majority of the objects in the plot is 
fully convective. While the lower activity limit in this plot is a detection limit, the upper 
limit is reliably determined and can be compared with published samples. In our sample, we 
obtain ~3x 10~ 4 , excluding the datapoint for TWA 10, which possibly is affected by a flare 
event (see §4.2p . For the mass range of our sample, this value is roughly comparable with 
the upper limi t in the Pleiades (IHodgkin et al.lll995l ). but clearly higher than in the Hyades 



1.4 x 10 4 , IStauffer et al.lll997al ). again indicating a decline of the general activity with 



age, as already discussed in §4.11 



As can be seen in Fig. 0, the upper limit of the range in activities is mostly flat. Thus, 
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activity is not strongly correlated with mass for 5kms _1 < vsini. This holds even when 
we only consider objects with radiative core and thus the potential to operate a solar-type, 
rotationally driven dynamo. It is also important to note that the distribution of rotational 
velocities for the non-active stars (not contained in Fig. [7J is indistinguishable from the active 
stars; they cover the full range from < 5 to 100 km s -1 , with a accumulation between 10 and 
20kms _1 . Moreover, among the four slowest rotators in our sample with vsini < 5kms _1 
(shown as upper limits in Fig. [7]), there is only one object with an activity level significantly 
below the range of datapoints for the faster rotators. Thus, by and large the rotation-activity 
correlation derived from Ha emission is flat in our sample. 

These resul ts can be compared phenom enolog;icalry with rot at ion- activity studies based 
on X-ray data. Ide la Reza fe Pinzonl ((200J) find that stars in TWA, BPMG, and TH are 
roughly comparable to T Tauri stars in the ONC in terms of their X-r ay properties. The ac 



tivity in the ONC h as been studied in detail in the COUP project fe.g. lFlaccomio et alj 
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Stassun et al.l 120041 ) . Both in the COUP data and in the sample of Ide la Reza fc Pinzon 
( 120041 ). there is no strong correlation between L x /L\> \ and rotation period. The rota- 
tion/activity relationship appears to be flat over a wide range of periods, interpreted as 
saturation with some indication for supersaturation, i.e. a decline of activity for the fastest 
rotators. This is very similar to what we observe in Ha. The two ultrafast rotators in Fig. [7J 
appear to have below average activity levels, which might be interpreted as supersaturation. 
The two additional ultrafast rotators not plotted in Fig. [TJhave no measurable activity level, 
thus confirming this trend. However, since we have only very few datapoints at high rota- 
tional velocities, this should be treated with caution. Still, it is interesting to note that the 
four ultrafast rotators are objects with radiative core, maybe implying that supersaturation 
might be associated to the presence of a solar-type dynamo. 

It is well established that field stars show a mostly linear relationship between rotation 
and relative X-ray luminosity (jRandichl 120001 ) . In young open clusters like IC2391, IC2602, 
and the Pleiades with ages ranging from 30 to 150 Myr, an intermediate situ ation is seen, with 
many objects in the saturated regime and an additional linear part (e.g. iPatten fc Simon 



19961 ). A hint of a linear relation might also be seen i n the sample of post T Tauri Lindroos 
stars with ages between 10 and 100 Myr analyzed by iHuelamo et al.l (120041 ) . Linear corre- 
lations bet ween X-ray flux and rotati o n rate have ad d itiona lly been found for young stars 
in Taurus (jStelzer fc Neuhauserl 120011 ). IStassun et al.l (120041 ) argued that the linear part of 
the rot at ion/ activity correlation in the ONC may be hidden in the objects for which no 
periods have been measured. However, studies of magnetic activity at very young ages are 
problematic, because accretion additionally affects both X-ray and Ha luminosities, which in 
principle requires the strict separation of accretors from non-accretors. Our Ha luminosity 
vs. v sini plot does not reveal a strong indication for a linear regime in the rotation/activity 
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relationship. The linear part cannot be hidden at low and thus unde tectable rotational 
velocities, as it has recently been found for field M stars ( Reinersl 2007 ). because even the 
slowest rotators in our sample show the same range of activity levels (with one exception, see 
above). In summary, it is still not clear if the linear part of the rotation/activity correlation 
is already established at ages < 30 Myr. 

In this context it is interesting to note that the rot at ion/ activity relation of young stars 
is similar to very low mass (VLM) objects with masses <O.3M . It is known that the ro- 
tational velocity at w hich saturation is reached drops quickly with decreasing object mass 



(Pizzolato et a 



(IDelfosse et al 



1120031). with the result that m ost VLM objects appear in the saturated reg: 
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lme 



Mohantv &: Basrill2003l ). As a consequence, the Skumanich type brak- 



ing law breaks down (Sills et al. 200ol ). and a weak exponential braking law is expected with 



f rot oc exp (— t) fjBarnesll2003bl ) . Such weak rotational braking is indeed req uired to model the 
rotat i onal e volution in the VLM regime on timescales of ~100-1 000 Myr (IScholz fc Eisloffel 
2004J . l2005h . 



In the canonical picture of the rotation/activity connection, the solar-type afl dynamo 
strongly depends on rotation, causing a linear relationship at low and moderate activity 
levels. The saturation effect is usually interpreted as an activity level where the stellar 
surface is covered by magnetic flux tubes and no further enhancement of activity (and rota- 
tional braking) is possible - thus the term 'saturation'. (It is unlikely that this corresponds 
to a surface completely covered with starspots, given the fact that many 'saturated' stars 
show strong photometric modulations due to rotation and thus have only a partially filled 
surface.) Both VLM objects and very young stars, however, are fully convective and thus 
cannot harbor a solar-type dynamo, which operates at the transition between convective and 
radiative zone (see JQ). The kind of alternative dynamo that produces t heir magnetic activ- 



ity, and how it depends on rotation, is still a matter of debate (see e.g. iDurney et al.lll993 



Chabrier fc Kukerll2006l ; iDonati et al.ll2006l ). But it is at least questionable to assume that 
the picture of the rotation/activity connection used for evolved solar-type stars can simply 
be extended to young stars and very low mass objects. For these types of objects 'saturation' 
can have two meanings: a) They are saturated and do not follow the Skumanich law, because 
they rotate too fast to be in the linear regime, as it is assumed in the standard paradigm, b) 
Their 'saturation' is the consequence of a magnetic field generation fundamentally different 
from solar-type stars. 

In the second case, 'saturation' is not merely the consequence of fast rotation, but a 
more fundamental sign of a change in the magnetic field generation (and thus the word 
'saturation' might be misleading). This idea has been proposed as an interpretation of 
rotation and activity data for open cluster stars by iBarnes J2003al jbll Basically, the Barnes 
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scheme suggest that the rotation/activity properties can be understood only in terms of 
the magnetic field generation: Fully convective obj ects in t he satu rated or supersaturated 
regime (on the 'C-sequence' in the nomenclature of iBarned (l2003al )) do not harbor a solar- 
type magnetic dynamo and thus do not follow the Skumanich type rotational braking law. 
As the stars evolve from the T Tauri phase to the ZAMS (and develop a radiative core), the 
fraction of objects on the C-sequence drops quickly an d reache s values < 10% at the age of 
the Hyades. While the quantitative predictions of the iBarnesI (j2003al ) scheme may not be 
convincing in all cases, the qualitative picture is consistent with the current rotation/activity 
data for the pre-main-sequence evolution. 

Many of our target stars in rj Cha, TWA, BPMG, and TH are too old to be fully 
convective (see Depending on their mass, they have developed radiative cores with 

substantial radii already. Thus, they present an interesting test case for magnetic field 
evolution. The fact that they show rotation/activity properties similar to younger stars 
(and to fully convective VLM objects) might indicate that it takes at least 30 Myr until 
the solar-type dynamo dominates the magnetic activity and rotational braking. This is 
supported by the weak rotational braking on the pre-main sequence found in §5.21 Future 
investigations of the magnetic field properties of pre-main-sequence stars as a function of 
age hold great potential to clarify these issues. 



6. Constraints on ages and radii: Rotation periods vs. vsiai 



By combining our measured v sin i with previously published rotational periods, we can 
derive stellar radii (times the unknown projection factor sini) as Rsini = (2it)~ 1 Pv sinz, 
and compare these to evolutionary models, to constrain the age. This gives an age estimate 
independent of other indicators such as the lithium abundance and color-magnitude diagram 
constraints, although all these estimates are dependent on the particular evolutionary model 
used. The estimated Rsini thus test the self-consistency of the models. 

Because of the unknown projection factor sini, a statistical sample is needed to derive 
the true R. Unfortunately, since our targets are widely distributed in the sky, monitoring 
campaigns to find photometric periods are time consuming, and only a small subset of our 
sample has photome tric periods measured. In total, 16 pe riods have been measured in r\ Cha 
( jLawson et al.l l200ll ) and TWA (jLawson fc Crausd 120051 ) . Of those, 13 are M dwarfs, out 
of which 12 have v sini above our detection threshold. In Fig. [8] we plot Rsini against 
effective temperatur e for those 12 targets, together with radius isochrones from models by 
(119981 ). The errors in the Rsini are entirely dominated by errors in v sini. 



Baraffe et al 
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Although the statistical sample (12) may seem small, we are helped by the fact that the 
probability density distribution for the projection fac tor sin % of a random ori entation favors 
close to edge-on geometries (see, e.g., Appendix A in lBrandeker et al.l (120061 )): 



/(sini) 



sm % 



sin 2 i 



Assuming a single age and the evolutionary models by iBaraffe et al.l (119981 ). equation 
(PQ) can be used to find a maximum-likelihood estimate for the age. To take into account 
the estimated measurement error a, and to mitigate the singularity at sini = 1, we assume 
the measured Rsini to be an outcome of a stochastic variable 1Z G -R mo d(£> T^)Y + E, 
where Y = sini is the projection factor distributed according to equation (CQ), E is normally 
distributed with zero mean and variance a 2 , and R mo d(t,T e s) is the model radius for a star 
of age t and effective temperature T e fj. The probability distribution of 1Z is obtained by 
numerical integration, 



fn(r\Rmod,^) = -5— I fy ( p— ) M r -x)dx 

1 f Rmod xexp[-(r - x) 2 /(2cr 2 )] 



RmodCrV^ir Jo J R? nd - x 2 



dx, (2) 



and the maximum likelihood by finding the maximum of the likelihood function 

C(t) = ^log f n (rj\R mod j, a j). (3) 
3 

To estimate conservative confidence intervals for this estimate, we integrate the proba- 
bility density function fn to get the cumulative probability function F-j^. We then find the 
age limits to and ^i such that the probabilities 



1 + a 

P(smi > ma,x{smij}\to) = — - — and (4) 

1 - a 

P(sini > max{sia ij}\ti) = , (5) 



where a is the significance and the probability function is 
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P(sinz > max{sini,}|t) = 1 - F n (R modtj max{sini i }|i? modii , <jj). (6) 

j 

Using the above relations we find the implied ages of r/Cha and TWA to be t^cha = 
13lg Myr and £twa = 9-1 Myr, respectively, where the quoted confidence interval is of 95% 
significance. These ages are slightly higher than, but consistent with, estimates from lit- 
erature (6 Myr for ^Cha and 8 Myr for TWA), indicating that the model radii with 95% 
confidence are good to within ~15%. 



7. Summary 

Rotation and activity are important parameters in the stellar pre-main-sequence evo- 
lution, because they trace changes of interior structure and magnetic fields as well as the 
dissipation of circumstellar disks. We present a spectroscopic study of rotation (measured 
as v sin i) and chromospheric activity (measured as Ha EW) for a sample of 74 young stars 
with spectral types F5-M5 in stellar associations with ages from 6 to 30 Myr. More than half 
of the objects are still fully convective, while the remaining fraction has already developed 
a radiative core. The analysis is based on an extensive set of multi-epoch high-resolution 
spectra obtained with the 6.5 m Clay Magellan telescope. We achieve a rotational velocity 
accuracy of < 5kms _1 . In the following, we summarize our results: 



The range and distribution of Ha EWs do not depend significantly on age in the 
considered age range; instead they are a strong function of spectral type (and thus 
stellar mass). Mid F to early K type stars have Ha in absorption, while most later type 
objects show emission. Until early K types, the Ha EW are mostly consistent with 
pure photospheric absorption, while for later spectral types chromospheric emission 
dominates. 

The spectral type at which Ha goes into emission in our sample is clearly earlier than 
in older clusters Pleiades and Hyades, but later than in very young T Tauri stars. This 
indicates a mass dependence i n the lifetime o f acti ve chromospheres. Using this as an 



age criterion, as suggested by lHawley et al.l (Il999l ). we find that the plausible age of 



TH is in the range between 10 and 40 Myr. 

The chromospheric activity measured in Ha clearly declines as a function of age from 
T Tauri stars (1-5 Myr) to post T Tauri stars in our sample (6-30 Myr) to ZAMS 
objects (50-100 Myr). 
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4. Many objects with spectral types later than early K show measurable variability in 
Ha EW on timescales of weeks and months, which can be attributed to chromospheric 
processes. 

5. Most objects in our sample have projected rotational velocities between 5 and 60 km s _1 . 
Additionally, four ultrafast rotators with vsini between 70 and 130kms _1 are seen, all 
in BPMG and TH. The maximum and minimum of the v sini range decreases between 
spectral types mid F to early K, indicating a dependence of rotation braking on the 
depth of the convection zone. 

6. The average rotational evolution between 5 and 30Myr is consistent with angular 
momentum conservation. It does not agree well with constant angular velocity i.e. 
' period-locking' . T his is the opposite of what has been observed for ages 1-5 Myr 



(iRebull et al.ll2004j ) and indicates a change in the rotational regulation at ages of ~5- 
10 Myr, coinciding with the average lifetime of (inner) disks. This may be interpreted 
with a scenario where the rotation is regulated by disk interaction at early ages, while 
they are free to spin up after the disks have disappeared. 

7. By comparing our data with rotational velocities in the Pleiades, we see some evidence 
for weak rotational braking on timescales of ~ 100 Myr. This might be an exponential 
or a Skumanich type rotational braking due to stellar winds (v oc t~ 1//2 ), which is 
switched on after the objects have developed radiative cores, i.e. after ~ 30 Myr. 

8. The rotation-activity relation, using Lua/Lhoh appears flat and thus 'saturated' in our 
sample. The maximum level of L-^ a /L^ i is ~ 3 x 10~ 4 , more or less independent of 
rotational velocity. There is no clear sign of a linear rotation/activity correlation at 
low vsini. 

9. The rotation- activity relation of stars with ages < 30 Myr is similar to fully-convective 
very low mass objects. The flat rotation/activity relation and the weak wind braking 
seen in these two object classes may not be due to 'saturation' of a solar-type rota- 
tionally driven dynamo, as suggested in the standard picture. Instead, the magnetic 
fields in these young objects are probably generated in a fundamentally different way 
from those in main-sequence stars. 

10. By comparing our rotational velocities with rotation periods from the literature, we 
find ages of 13^g Myr and 9lJ Myr for rj Cha and TWA, respectively, consistent with 
previous estimates from other methods. Th i s agre ement indicates that the stellar radii 



for M dwarfs from models by lBaraffe et al.l (119981 ) are good within ~15%. 
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Table 1. Summary of results for r] Cha. 



Star 


Ha EW 


Ha EW a 


v sin i 


vs'mi a 


Spectral Type a 


rj Cha 3 


-1.99 


0.17 


10.50 


0.67 


M3.25 b 


r?Cha 4 


-3.40 


0.61 


5.96 


1.00 


K7 a 


rj Cha 5 


-8.57 


4.29 


8.75 


1.48 


M4 b 


r\ Cha 6 


-5.04 


0.42 


20.89 


1.05 


M2 a 


•q Cha 10 


-1.15 


0.28 


<5.0 


0.93 


K7 a 



Note. — The name of the star, Ho equivalent width, standard deviation in 
the Ha equivalent width, v sin i, and standard deviation in v sin i are provided. 
A positive Ha EW denotes absorption. 

HZuckerman fc Sond j2004) 

iLuhman fc Steeehsl ll2004h 



Table 2. Summary of results for TWA. 



Star 


Ha EW 


Ha EW 


a 


v sin i 


v sin i 


<7 


Spectral Type 


TWA 2a 


-1.84 


0. 


21 


12.78 


1. 


.00 


M2e b 


TWA 3b 


-6.14 


0. 


.89 


12.20 


6. 


.59 


M3.5 b 


TWA 7 


-5.82 


0. 


.85 


<5.0 


1. 


.59 


Ml a 


TWA 8a 


-8.00 


1. 


42 


<5.0 


1. 


.15 


M2 b 


TWA 8b 


-13.3 


1. 


84 


11.20 


2. 


.52 


~M5 C 


TWA 9a 


-2.14 


0. 


.45 


11.26 


0. 


.53 


K5 b 


TWA 9b 


-4.31 


0. 


60 


8.39 


0. 


.61 


Ml b 


TWA 10 


-13.6 


9. 


63 


6.33 


1. 


.18 


M2.5 a 


TWA lib 


-3.45 


0. 


62 


12.11 





.93 


M2.5 b 


TWA 12 


-4.83 


0. 


85 


16.22 





.94 


M2 b 


TWA 13a 


-3.01 


0. 


69 


10.46 


1 


.10 


Mle b 


TWA 13b 


-3.00 


0. 


74 


10.33 


1. 


.19 


M2e b 


TWA 15a 


-8.81 


0. 


51 


21.33 


2 


.16 


M1.15 b 


TWA 15b 


-8.59 


1. 


35 


32.33 


1. 


.86 


M2 b 


TWA 18 


-3.32 


0. 


36 


24.12 


0. 


.74 


M0.5 b 


TWA 19a 


2.60 


0. 


13 


28.35 


0. 


.71 


G5 b 


TWA 21 


-0.11 


0. 


18 


6.00 


1. 


.00 


Ml b 


TWA 22 


-11.5 


1. 


.95 


9.67 


2 


.37 


M5 a 


TWA 23 


-2.44 


0. 


19 


14.78 


1. 


.72 


Ml a 


TWA 24s 


-0.30 


0. 


16 


13.87 


1. 


.13 


K3 a 


TWA 25 


-2.36 


0. 


54 


11.78 


1. 


.86 


M0 a 



Note. — See notes for Table [T] for column heading explanation. 
' IZuckerman fc Sond (120041 1 
y<le la Reza fc Pinzonl 1120041 ) 

c Spectral type approximated by visual comparison with other stars in this 
group. 
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Table 3. Summary of results for BPMG. 



Star Ha EW Ho EW a v sin i v sin i a Spectral Type a 



AO Men 


-0.56 


0.12 


15.96 


0.69 


K6/7 


AU Mic 


-2.24 


0.81 


8.49 


0.97 


Ml 


HD 164249 


5.43 


0.09 


22.50 


1.97 


F5V 


HD 181327 


5.09 


0.05 


20.83 


1.33 


F5.5 


HD 35850 


4.31 


0.18 


52.00 


1.00 


F7 


PZ Tel 


2.29 


0.17 


77.50 


2.81 


KOVp 


V343 Nor 


2.22 


0.15 


17.00 


1.00 


K0V 


GJ 3305 


-2.15 


0.27 


5.30 


1.03 


M0.5 


GJ 799n 


-10.9 


1.61 


10.56 


2.13 


M4.5e 


GJ 799s 


-9.25 


0.61 


17.00 


3.54 


M4.5e 


HIP 23309 


-0.77 


0.13 


5.77 


0.73 


M0.5 


HIP 23418A 


-6.55 


0.00 


7.67 


2.08 


M3V 


HIP 23418B 


-6.06 


0.00 


21.00 


4.36 


~M3 b 


HIP 112312 


-6.61 


0.14 


14.00 


1.73 


M4e 


HIP 112312B 


-8.16 


0.46 


24.33 


4.93 


M4e 


HIP 12545 


-0.56 


0.00 


9.30 


0.64 


M0 



Note. — See notes for Table [T] for column heading explanation. 
HZuckerman fc Sond J2004h 

b Spectral type approximated by visual comparison with other stars in this 
group. 
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Table 4. Summary of results for Tuc-Hor 



otar 


net rLW 




v sin i 


v sin i a 


Spectral Type a 


CD -53544 


-1.43 


0.24 


82.22 


5.02 


K6Ve 


CD -60416 


-0.48 


0.04 


10.11 


0.60 


K3/4 


CPD -64120 


-0.19 


0.12 


30.22 


1.09 


KlVe 


CSC 8056-0482 


-5.31 


0.44 


34.22 


5.61 


M3Ve 


CSC 8491-1194 


-4.14 


0.36 


12.78 


0.83 


M3Ve 


CSC 8497-0995 


-0.63 


0.22 


6.56 


0.53 


K6Ve 


HIP 107345 


-1.41 


0.15 


6.44 


1.33 


Ml 


HIP 1993 


-1.01 


0.11 


9.50 


0.55 


Ml 


HIP 2729 


-0.68 


0.26 


127.5 


3.94 


K5V 


HIP 3556 


-0.79 


0.02 


< 5.0 


0.52 


M3 


HD 13183 


2.70 


0.13 


21.00 


0.71 


G5V 


HD 13246 


4.13 


0.13 


29.78 


1.79 


F8V 


HD 8558 


2.62 


0.10 


12.11 


0.60 


G6V 


HD 9054 


0.76 


0.11 


< 5.0 


0.53 


K1V 


HIP 105388 


2.58 


0.10 


12.83 


0.39 


G5V 


HIP 107947 


4.26 


0.05 


30.42 


1.38 


F6V 


HIP 108422 


1.55 


0.54 


139.8 


9.42 


G8V 


HIP 1113 


2.57 


0.05 


6.00 




G6V 


HIP 1481 


4.06 


0.00 


22.67 


0.52 


F8/G0 


HIP 16853 


3.41 


0.10 


17.80 


0.41 


G2V 


HIP 21632 


3.02 


0.00 


17.53 


0.64 


G3V 


HIP 22295 


3.54 


0.29 


41.33 


2.10 


F7V 


HIP 30030 


3.31 


0.17 


40.40 


1.24 


GO 


HIP 30034 


1.30 


0.26 


9.73 


1.71 


K2V 


HIP 32235 


2.72 


0.05 


10.25 


0.45 


G6V 


HIP 33737 


0.77 


0.52 


8.93 


0.26 


K3V 


HIP 490 


3.72 


0.03 


14.50 


0.55 


GOV 


HIP 9141 


3.03 


0.07 


14.78 


0.44 


G3/5V 


TYC 5882-1169 


0.74 


0.16 


6.80 


0.56 


K3/4 


TYC 7065-0n 


1.45 


0.02 


22.33 


1.03 


K4V 


TYC 7065-0s 


1.58 


0.06 


14.00 


0.89 


~K4 b 


TYC 7600-0 


1.76 


0.25 


18.80 


0.86 


Kl 



Note. — See notes for Table [T] for column heading explanation. 
HZuckerman fc Sond j2004h 

b Spectral type approximated by visual comparison with other stars in this group. 
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Fig. 1. — Standard deviation of v sinz versus average vsini. Only 4 objects out of 74 (5%) 
have standard deviations greater than 5kms _1 (see §3.31) . 
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Fig. 2. — Absolute differen ce between measured usinz and literature values compiled by 
de la Reza fc Pinzonl ((200J). The error bars correspond to the scatter in our multi-epoch 
data. With the exception of HIP 2729 and PZ Tel, the deviations are not larger than 8 km s -1 
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Fig. 3. — Ha EW as a function of spectral type. FO corresponds to 30, MO to 60. The dotted 
line shows the Ha EW for non-active field stars; the dashed line marks the zero level. 
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Fig. 4. — Ha EW standard deviation as a function of spectral type compared with the 
measurement error (dashed line). FO corresponds to 30, MO to 60. The datapoint for 
TWA 10 at spectral type M2.5 and Ha EW a of 9.6 A is not plotted. 
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Fig. 5. — Projected rotational velocities v sini as a function of spectral type. FO corresponds 
to 30, MO to 60. The four objects above the dashed line are the so-called ultrafast rotators 
in our sample. 
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Fig. 6. — Projected rotational velocities vsini as a function of age, symbols as in the pre- 
vious figures. Upper panel: Big oct agons mark the med ian values in the four associations 
and the average at 5 Myr (taken from iRebull et al.l (12004 )). The solid line shows the evolu- 
tion assuming constant angular momentum (i.e. spin-up, model A), the dashed line constant 
angular velocity (i.e. period locked, model B). Lower panel: In addition, the median vsini 
in the Pleiades is plott ed (big octag o n at 125 Myr), together with the quartiles (horizon- 
tal lines), derived from iQueloz et al.l (119981 ). The solid lines show the rotational evolution 
assuming constant angular momentum (as in the upper panel, model A), the dotted lines 
assume angular momentum losses following a Skumanich type law (v oc t~ 1//2 , model C). 
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